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1. Introduction 

The idea that the strong interactions are described by a string theory which is in 
some sense dual to perturbative QCD is a major challenge for high energy theory. More 
generally, a D-dimensional confining Yang-Mills theory is expected to define a string theory 
with D-dimensional target space stable in the interval 2 < D < 4. 

This "YM string" has been constructed so far only for D = 2 |JJ,@. It has been 
shown that the partition function of pure U(N) gauge theory defined on a two-dimensional 
manifold of given genus and area can be represented in terms of a weighted sum over maps 
from a worldsheet E^ to spacetime E^. The allowed worldsheet configurations represent 
minimal area maps Y*w — ► Ey. The latter condition is equivalent to the condition that the 
embedded surfaces are not allowed to have folds. Ramification points are however allowed, 
as it has been suggested in earlier studies ||. In other words, the path integral of the 
' YM2 string" is over all branched covers of the target manifold 0,0. 

Unfortunately, this construction is highly involved and does not easily reduce to a 
system of simple geometrical principles^. One can nevertheless speculate, using the analogy 
with the random- walk representation of the O(N) model, that there exists an underlying 
2D string theory with clear geometrical interpretation. The YM2 string is obtained from 
the latter by tuning the interactions due to ramification points and adding new contact 
interactions via microscopic tubes, etc. 

One is thus led to look for the most general 2D string theory whose path integral is 
given by branched covers of the target space. Any such string theory is invariant under 
area-preserving diffeomorphisms of the target space Ey. Its partition function depends on 
Tit only through its genus and area. In addition to the topological coupling constant iV -1 
there is a set of couplings t n associated with the interactions due to ramification points of 
order n. 

In this letter we propose a discretization of the path integral for such string theories, 
in which the target space Et represents a two-dimensional simplicial complex with given 
topology. The branch points (the images of the ramification points) are thus located at 
the vertices of the target lattice Ey. We will not try to establish a worldsheet description 
but instead construct and solve a class of equivalent matrix models. Our approach is 

1 The random surface representation of YM2 on a lattice found in [[| is geometrically clear 
but has the inconvenience of being highly redundant. For example, the folds are not forbidden 
but their contribution vanishes as a result of cancellations. 
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thus analogous to quantizing Polyakov string theory by employing randomly triangulated 
surfaces. 

The matrix model associated with the target space £t will be formulated as a 2D 
lattice gauge theory whose local link variables are complex N x N matrices. However, 
the matrices are no longer unitary: the Haar measure on U(N) is replaced by a Gaussian 
measure. The model resembles very much the 2D Weingarten model with the important 
difference that the lattice action represents a sum over the positively oriented cells only. 
The latter restriction eliminates from the string path integral all surfaces containing folds 
(which are believed to cause the trivial critical behavior of the standard Weingarten model 
0). We will demonstrate that these 2D matrix models are exactly solvable for any target 
space lattice and any N by employing the character expansion methods recently developed 



i n 5 B 5 5 E3 • We pay special attention to the cases of spherical and toroidal topology. 
In the first case we observe that the string theory exists only for a target space with 
sufficiently large area. At the critical area a third-order transition takes place due to the 
entropy of the ramification points. In the second case we find the same partition function 



as the chirally perturbed conformal field theory considered recently by R. Dijkgraaf [11 



In this letter we restrict our attention to the discrete case since it involves combinatorial 
problems interesting on their own. The continuum limit of an infinitely dense target lattice 
will be considered in detail elsewhere [O . 



2. Definition of the models 

By target space we will understand an oriented triangulated surface (two- 
dimensional simplicial complex) of genus G containing Ao points, J\T\ links and M2 two- 
dimensional cells. These numbers are related by the Euler formula 

AT - M + M 2 = 2 - 2G. (2.1) 

In addition, each cell c is characterized by its area A c . 

We will first consider the simplest string theory with this target space, in which all 
ramification points have Boltzmann weight one, and the corresponding matrix model. The 
generalized model, to be considered in sect. 6, will depend on a set of external field variables 
associated with the points of Sy. 
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At each link I is denned a field variable $^ representing an iV x N matrix with complex 
elements. The partition function of the matrix model is defined as 

z = f n ex p(^ n Tr$c )' ( 2 - 2 ) 

J i c 

where (3 C = e~ Ac , $ c denotes the ordered product n^edc ^ °f variables along the 
oriented boundary dc of the cell c, and the integration over the link variables is performed 
with the Gaussian measure 

[V* t ] = {NM N *ll ij d* ij d** iJ e- N ^ *<**. (2.3) 

Note that we omit the complex conjugate from the plaquette action, eliminating thereby 
orientation reversing plaquettes. This is why we call the model with partition function ( |2.2j ) 
a chiral matrix model. The model is invariant under complex conjugation of the matrix 
variables and reversing the orientation of the target space. 

The perturbative expansion of ( |2.2[ ) results in a representation of the free energy 
T = InZ in terms of connected lattice surfaces Hw embedded in the target surface Ey. 
Each of these surfaces is obtained by placing plaquettes on the faces of the lattice Ey and 
gluing any two together along the edges. All plaquettes should have the same orientation, 
which means that the surfaces cannot have folds. It is geometrically evident that the 
number of plaquettes, say n, covering each cell is constant throughout Ey. Thus the 
surface Evk is wrapping n times E j- and its area is A\y = uAt where At = J2 C A c is the 
total area of the target space. The surface E^ may have ramification points whose images 
are points of Ey. The map Evk — > Ey defines at each point p G Ey a branching number 
B p . The branch points are the points with B p ^ 0. By the Riemann-Hurwitz formula 
2(7 — 2 = n(2G — 2) + ^2 p B p where g is the genus of the surface Evk- The free energy 
JF = InZ of our chiral model can be written as 

oo oo 

F=Y1 e ~ nAT n2 ~ 29 F(^ 9\G, ATo) (2.4) 

n=l g=G 

where J-(n, g\G, Ao) is the number of no-fold surfaces of genus g covering n times Ey. This 
number depends on Ey only through its genus G and the number of the A/o of allowed 
locations of the branched points. 
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3. Exact solution by the character expansion method 



Applying to ( |2.2j ) the same strategy as in ref. |I(J, we expand the exponential of the 
action for each cell c as a sum over the characters x h °f the polynomial representations 
of GL(N). We use the shifted weights h — {hi, h 2 , . . . , h^} where hi are related to the 
lengths mi, mjv of the rows of the Young tableau by hi = N — i + rrii and are therefore 
subjected to the constraint hi > h 2 > . . . > h^ > 0. We will denote by \h\ = T>irrii the 
total number of boxes of the Young tableau and by 

the dimension of the representation h. The character expansion of the exponential then 
reads 

h 

Now the model can be diagonalized in /i-space by using the following two simple but 
powerful facts. Firstly, one easily proves J7|,PD| 



where we have defined 



[V$] Xh (^)=n h A h , (3.4) 

where flh is given as above by ( |3.3| ). Secondly, one uses that complex matrices may be 
diagonalized by a bi-unitary transformation U<&dV ] ' , where U, V are unitary and $d is a 
positive definite and diagonal matrix. Then one uses the formula ( |3.4[ ) together with the 
fission and fusion rules for unitary matrices 

/ VU xMUBW) = J- xM) X h (B), 
J VU Xh (AU) Xh , (WB) = 8 h>h , i- Xh (AB), 
to derive fission and fusion rules for complex matrices corresponding to the measure ( |2.3|) : 

f[VQ] Xk (A*B&) = ^ Xh {A) Xh (B), 

r n (3 ' 6) 

J [V*\ X h (A®) X,A&B) = 8 h , h , -± Xh (AB). 



Therefore, complex matrices behave exactly like unitary matrices apart from the extra 
local factor for each integration!. 

Having at hand (|3.6|), we can proceed with the exact solution of ( |2.2|) , just as in the 
gauge theory case |13j , |14]] . Using the fusion rule one progressively eliminates links between 
adjoining cells until one is left with a single cell. The remaining links along that plaquette 
are integrated out by applying the fission rule. One arrives at the formula 

Z = Y, A l~ 2G ^h° +2G ~ 2 e~ WAT . (3.7) 

h 

We observe that, in accord with the surface representation ( |2.4|) , the solution's only de- 
pendence on the original target space lattice data is through the genus G, the number of 
vertices Ao, and the total area At of the target space. It is therefore clear that a given 
character expansion ( |3.7| ) corresponds to many different complex multi matrix models ( |2.2| ) 
with the same partition function: As many, as we can form simplicial complexes differing 
in their internal connectivity but constrained by the global data G, Ao and At- (This 
is a lattice version of the invariance with respect to area- preserving automorphisms.) In 
particular we can represent the target space by a one-cell complex and write down a matrix 
model with a minimal number Ai = Ao — 1 + 2G of link variables giving the expansion 



„M)-1+2G / G N -l+2G \ 

Z= n [V* e ] exp le~ AT iV Tr [ J] $ 2s -i^L-i$L II M 3 " 8 ) 

* s=l fc=2G+l ' 



Note that the trace in the exponent of (|3.8|) corresponds to the element of the ho- 
motopy group of a surface with G handles and Nq punctures which is equivalent to the 
identity (Fig. 1). 



2 The fission and fusion rules actually hold for any U(N) invariant measure if we define 

the factor Qh by the integral ( |3.4D . In particular one could consider potentials higher than linear 



in in (|^3|) . 
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Fig.l: The target space of the effective one-plaquette model. 



4. Spherical target space 

In the case of a spherical target space (G = 0) the partition function ( |3.7| ) becomes: 

Z G=Q = J2^l Of" 2 e-l^. (4.1) 
h 

The simplest example of Ao = 2 is described by a Gaussian integral. There is only one 
surface covering n times a sphere with two punctures and it contributes \&~ nAT to the 
free energy (the factor — coming from the cyclic symmetry). No higher-genus surfaces are 
contributing. Summing on n gives therefore 

T = -N 2 \n(l-e~ AT ), (4.2) 

which is also the result of the Gaussian integration. 

If we have Mq > 3 the situation becomes non-trivial. Now there are contributions from 
non-spherical (g > 0) world sheets as well. Concentrating on the leading (g = 0, i.e. N = 
oo) term, we can proceed by using the saddlepoint techniques of 0. Here the sum (fOj) 
describes a "gas" of mutually repulsive weights at thermal equilibrium and is dominated by 
the contribution of a "classical" configuration {h}. One introduces continuum variables h = 
jfjr and a density p(h) = -h ^ S(h — hi). The saddlepoint density is found by functionally 
varying (|4.1| ) with respect to h; this gives 



It is straightforward to explicitly solve this equation and calculate the free energy. Here 
we will merely write out the equation determining the critical behavior; one finds, for the 
variable u = \{\fa, + Vb) 2 , the algebraic equation 

(JVo-1) e~ AT u^- 1 =u-l. (4.4) 

This equation predicts for J\f > 3 a critical area At = A T > 0: The sum over surfaces 
diverges if the area of the target manifold falls below A T . A quite similar transition has 
already observed in YMi by Douglas and Kazakov [IS]. One finds for the free energy near 
the critical area 

T ~ N 2 (A T - A c T )i (4.5) 

i.e., the same critical behavior as the Polyakov string in zero dimensions. We can under- 
stand the critical behavior by considering e.g. Ao = 3 

Z = J [V$i][V$ 2 ] e qN i**! (4.6) 

where q = e~ Ar . Eq. (|4.6|) generates planar diagrams which can be interpreted as abstract 
plaquettes without embedding. Its critical behavior is in the university class of pure gravity. 

For Ao > 4 the models ( ft.8|) are intractable without our character expansion tech- 
niques. The qualitative behavior remains the same: As seen from (|4.4|) , we can always find 
the behavior of pure 2D gravity. The critical size of the target sphere is given by 

'M> - 1' 
JV -2. 

and tends to infinity in the continuum limit Ao — >• 00. It is therefore evident that the 
theory should be modified in order to have a sensible continuum limit. 



A T = 2 KM, - 1) + (A- - 2) In (^l) ■ (4-7) 



5. Toroidal target space 

If the target space is a torus (G = 1) the Vandermonde determinants in ( |3.7|) cancel, 
the free energy is now order O(N ) and the saddlepoint method is no longer applicable. 
However, we can do even better in this case and immediately give a result to all orders 
in -^j. Each term in the sum ( |3.7| ) factorizes into a product over the rows of the Young 
tableau; the sum remains nevertheless non-trivial due to the ordering constraint on the 
weights. Using ( |3.3|) , rewrite Z as 

2ci=e n (i-W- (5 - x) 

h {i,j}eh 



7 



Here denotes the a box in row % and column j and the product goes over all boxes in 

the tableau h. Slicing the tableau through the diagonal and counting the fraction of boxes 
in the rows of the upper half and the columns of the lower half we can elegantly express 



(0) as 



Z, 



G=l 



dz 
2niz 



oo n 7 Kj 



n=0 



k=l 



1 q n+ 2 



n 

k=i 



k \M) 

n) J 



(5.2) 



where we have denoted q = e~ Ar . The contour integral ensures that an equal number of 
rows and columns emanate from the diagonal; it eliminates configurations that cannot be 
interpreted as a tableau!. The leading order of the free energy gives the partition function 
of the noninteracting string. It is immediately extracted from either ( |5.1|) or (|5.2|) and is 
identical to what one finds for the continuum Yang-Mills theory on a torus 



T = - ln (! - 9 n ) + 0{N~ 2 ). 



(5.3) 



n=l 



The differences start with the terms. These are in principle computable in terms of 

IT 



MM 



quasi modular forms as in |T 

Similarly to the YM2 theory, the matrix models with toroidal target spaces can be 
interpreted as chiral deformations of a two-dimensional free field theory defined by the 
action CT,p^H,pl[ 



S 



n + 1 



(dip) 



n+l 



(5.4) 



The partition function (|5.2[) is equal to the zero U(l) charge sector of the partition function 
of (|5.4|) computed on a torus, with an appropriate choice of the couplings {s n }. 

Let us mention an interesting interpretation of the model with A/o = 1 whose target 
space is a torus with one puncture. Explicitly, we have 



Z = J [D$ 2 ] e qN i**i* a *I< 



(5.5) 



The model is superficially very similar to (|4.6|) of the last section. Again a network of 
square plaquettes (corresponding to the vertices of ( |5.5|) ) is generated. But here a local 



3 It is interesting to point out that a very similar result was first found by Douglas [16] in the 
context of continuum 2D Yang-Mills theory on the torus. The combinatorial derivation of (|5.2|) 



is due to Dijkgraaf [17]. 
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rule suppresses all positive curvature: The coordination numbers at the vertices of such a 
surface can take values 4, 8, 12, which correspond to zero or negative curvature. No local 
curvature fluctuations, leading to pure quantum gravity behavior, are possible: For fixed 
genus g the surfaces contain only a finite number of quantized, negative curvature defects 
of total curvature — 47r(g — 1). The model is therefore similar in spirit to the "almost flat 
planar graphs" of f|], where discrete surfaces of vanishing or positive local curvature were 
considered. 



6. The general case 

Let us now consider the discretized string theory in the general case when a ramifica- 
tion point of order n = 2, 3, 4, ... associated with the point p £ Ej- is weighted by a factor 

(p) 

t n ■ It is straightforward to appropriately modify the original matrix theory and solve it 

by repeating the same steps explained in section 3. 

(v) 

We will parametrize the couplings tm by an N x N external matrix field 

t« = j^Tr(B p r (6.1) 

and insert the matrix B p into the r.h.s. of ( |2.2|) . This should be done in the following way: 
For each point p we choose one of the cells c such that p £ dc and insert the matrix B p 
into the corresponding trace. For example, the one-plaquette theory (|3.8|) becomes 



„M)-1+2G / G M)-1+2G 

Z= J] [V* e ] exp e~ A - N Tr[i? M) J] * 2 ._ 1 * 2 .*L-i*L II B ^^l 

J 1=1 V s=l fc=2G+l 

(6.2) 

Then one can proceed as in sect. 3, arriving at the following final expression replacing 

The model ( |6.3| ) can be solved in the spherical limit using the methods developed in 
0-Q . Let us note that in the case of a sphere with three punctures (G = 0,Ao = 3) one 
reproduces the model of the dually weighted planar graphs in the most general formulation 



presented in [10 



2(a t =o,g=o,at =3) = X 1 X h (Bi)x h {B2)x h (B 3 ). (6.4) 

h h 
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7. Concluding remarks 

1. We have constructed a discretization of the most general string field theory of 
nonfolding surfaces immersed in a two-dimensional compact spacetime. In the continuum 
limit the sum over the positions of a branch point should be replaced by an integral with 
respect to the area. This means that the couplings t n should scale as 

where t n are the interaction constants of the continuum theory. The properties of the 



expansion (|6.3| ) in this limit will be studied elsewhere JT! 

2. The YM-string should be obtained by a special tuning of the coupling constants and 
by introducing contact (tube-like) interactions. These interactions can be implemented by 
considering the 5-matrices as dynamical fields. 

3. We point out an interesting equivalence between ensembles of coverings of a sphere 
with fixed number of punctures and an ensemble of abstract (i.e., without embedding) 
planar graphs. It is therefore not a miracle that we find a third order phase transition as 
in the case of pure 2D gravity. In our case the transition is due to the entropy of surfaces 
with a large number of ramification points. A very similar transition has already observed 



in YM2 by Douglas and Kazakov [15 



4. The chiral matrix models defined on a torus give a nice geometrical interpretation 



of the general chiral deformation of topological theory studied by R. Dijkgraaf fil]] . The 
sum over surfaces can again be interpreted in terms of abstract planar graphs describing 
random surfaces with non-positive local curvature. 
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